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ABSTRACT 

The paper presents a theoretical formulation based on polynomial shape function and application of Ritz method. 
In this study, the free vibration of simply supported panel with one free edge was analyzed. The Polynomial shape function 
derived was substituted into the potential energy functional, which was minimized to obtain the fundamental natural 
frequency. This research focused on aspect ratio from 0.1 to 2.0 with 0.1 increments. The values of fundamental natural 
frequencies of the first mode obtained were compared with those of previous research works. The present values for aspect 
ratio 0.5, 1.0 and 1.3 were 4.9469, 12.3435 and 19.1582 respectively. 
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1. INTRODUCTION 

There has not been any exact solution on the vibration of rectangular panels with one free edge. Over the years, 
problems have been treated by the use of trigonometric series as the shape function of the deformed panel or by using 
method of superposition. Research has been carried out on the problems from equilibrium approach and others solved the 
problems from energy and numerical approaches. However, no matter the approach used, trigonometric series has been the 
most widely used shape functions. The analysis of the free vibration of panels was well documented by Leissa (1973) and 
includes a variety of boundary conditions and aspect ratio using trigonometric series. Gorman (1982) solved problem on 
free vibration of rectangular panel with different boundary conditions, aspect ratio and Poisson ratios using method of 
superposition. 

2. FORMULATION OF FUNDAMENTAL NATURAL FREQUENCY 

Chakraverty (2009) gave the maximum kinetic energy functional as 

K max = ^JJphW 2 {x,y-)dxdy (1) 
Making use of the non dimensional parameters, R and Q, equation (1) becomes 

= a -^H W i d RdQ (2) 

Where p is the weight per unit area of the plate, h is the plate thickness and X is the frequency of panel vibration. 
The maximum strain energy functional for a thin rectangular isotropic panel under vibration was given by Ibearugbulem 
(2012) as follow: 
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as: 



Umax = j/J l(W" x ) 2 + 2{W" x y) 2 + (W"y) 2 ]dxdy (3) 
Adding equations (2) and (3) gave the total potential energy functional of rectangular panel under lateral vibration 



Umax = ^ JJ [±(W R ) 2 +^- 2 {W" m ) 2 + ^(W" Q ) 2 ] dRdQ - ^ jf w 2 dRdQ 
Factorizing out b/a 3 gave: 

Umax =^/j[(W' B ) 2 + 2 -£(W" R iy +£(w"Qy] dRdQ _^lj jw 2 dRdQ 
IfP= a/b, then: 

Ilmai=J // [(W" fl ) 2 + 2P 2 (W" R( iy + P 4 (W"Q) 2 ]d/?dQ - ^±fjw 2 dRdQ 
If p= b/a, then: 

Umax=^fi[(W" R ) 2 +^(W'^) 2 +^(W'^ 2 ]dRdQ-^Jjw 2 dRdQ 
From equation (4), if a/b 3 is factorized out then: 

Umax = ^jj g(W"«) 2 + 2 -f(W"^) 2 + (W'Q) 2 ] dRdQ - ^ jj w 2 dRdQ 
IfP = a/b, then: 

nmax = ^_tf{l_ iw >>Ry +^(W'W) 2 + ( W "Q) 2 ]dRdQ- a -^tfw 2 dRdQ 
If P = b/a, then: 

Umax = g^J [P 4 (W" fl ) 2 + 2P 2 (W" fl Q) 2 + (W"<2) 2 ]dRd<2 - ^±fJw 2 dRdQ 

T A YLOR-MCL AURIN ' S SERIES SHAPE FUNCTION 
Displacement Function for SSFS and SSSF Panel 
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Figure 1 

Deflection Equation 

The deflection equation, according to Ibearugbulem (2012), is: 
W = (a„ + a^ + a 2 R 2 + a 3 R 3 + a 4 R 4 )(b + biQ + b 2 Q 2 + b 3 Q 3 + b 4 Q 4 ) 



(11) 
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Meaning that W = W R . W Q (1 2) 

Boundary Condition (for panel A) Boundary condition (for panel B) 

W(Q = 0) = 0, W(R = 0) = W(R = 0) = 0, W(Q = 0) = 

W"(Q = 0) = 0,W"(R = 0) = W"(R= 0) = 0,W"(Q = 0) = 

VR(Q = 1) = 0, W(R = 1) = VR(R = 1) = , W"(Q = 1) = 

For simply supported edges, both kinematic and dynamic boundary conditions were satisfied. At the free edge, the 
kinematics boundary conditions were satisfied. The dynamic boundary conditions are the shear force and bending moment. 
Only the shear force condition was satisfied at the free edge. Using panel B, these gave: 

W R =a 4 (8R-4R 3 +R 4 ) (13) 

For W Q 

W Q = b 4 (Q-2Q 3 +Q 4 ) (14) 
Remembering from that 
w = W R . W Q 

W = A (8R-4R 3 +R 4 ) (Q-2Q 3 +Q 4 ) (15) 
The above equation is the displacement function for SSSF panel while the displacement equation/function for 
SSFS was W = A (R-2R 3 +R 4 ) (8Q-4Q 3 +Q 4 ) (16) 
Total Potential Energy for SSSF Panel 

Integrating the squares of the differential equation (15) with respect to R and Q 



J j (W" R ) 2 dRdQ = 4 2 (76.8)(0.04920634921) = 3.779047619 A 2 

o o 

l l 

J j (W' R ' Q ) 2 dRdQ = ,4 2 (31.08 5 71429)(0.48 5 7148 5 7) = 15.09 8 7755 L4 2 

o o 

i i 

j j (W" Q ) 2 dRdQ = 4 2 (12.5968254)(4.8) = 60.464761924 2 



o o 
l l 



| j (W) 2 dRdQ =4 2 (12.5968254)(0.04920634921) = 0.6184378964 2 

o o 

Substituting the results into integral (6), (7), (9) and (10) 
IfP = a/b 

Umax = ^[3.779047619 + 30.19755102P 2 + 60.46476192P 4 ] -^^[0.6198437896] (17) 
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P= b/a 

n max = [3.779047619 + 3 J^£1 + ™ ^™A _ [0.6198437896] (18) 

2s L P P J 2 

IfP = a/b 

„ DaA 2 [3.779047619 30.19755102 ,„-,-_,.. „„1 ab^ 2 phA 2 r „ , . . _ _„„ , n ,„ 

11 ma^ = —j- 1 — + — + 60.46476192 1 — [0.6198437896] (19) 

IfP=b/a 

I! max = ^0- [3.779047619P 4 + 30.19755102P 2 + 60.46476192] - ab ^ 2phA2 [0.6198437896] (20) 

Minimizing (17), (18), (19), (20) by - ^g™"* = ® anc ^ ma king ^ subject of formula. Table 1, 2, 3 and 4 gives the 
results of fundamental natural frequency and the comparison of present result with other research works. 

Figure 1 to 8 shows the graphical model of parameter (k) of fundamental natural frequency results, for both 
X-X and Y-Y axis 
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Figure 1: Graph of SSSF Panel, P=a/b (Note: y = k and x = Aspect Ratio, p) 

The free edge of the panel is on X - X axis. With respect to length (a), the natural frequency increases as the 
aspect ratio increases, the polynomial equation curve is represented by y =9.878x 2 - 0.004x + 2.469. 
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Figure 2: Graph of SSSF Panel, P=b/a (Note: y = k and x = Aspect Ratio, p) 
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The free edge of the panel is on X - X axis. With respect to length (a), the natural frequency decreases as the 
aspect ratio increases, the power equation curve is represented by y = 13.68x"' 75 




Figure 3: Graph of SSSF Panel, P = a/b (Note: Y = k and x = Aspect Ratio, p) 

The free edge of the panel is on X - X axis. With respect to width (b), the natural frequency decreases as the 
aspect ratio increases. This model consist of one polynomial equation y = 23436x 6 - 87611x 5 + 13164x 4 - 10152x 3 + 
42393x 2 - 9161.x + 837.9 which stop at P= 1 and linear equation y = -1.505x + 13.35. 
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Figure 4: Graph of SSSF Panel, P = b/a (Note: y = k and x = Aspect Ratio, p) 

The free edge of the panel is on X - X axis. With respect to width (b), the natural frequency increases as the 
aspect ratio increases, the polynomial equation curve is represented by y =0.001x 3 + 2.465x 2 + O.OOOx + 9.876. 
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Figure 5: Graph of SSFS, P = a/b (Note: y = k and x = Aspect Ratio, p) 
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The free edge of the panel is on Y - Y axis. With respect to length (a), the natural frequency increases as the 
aspect ratio increases, the polynomial equation curve is represented by y = O.OOlx 3 + 2.465x 2 + O.OOOx + 9.876. 



- 23436X 6 - S7611X 5 4- 13164**- 10152k 3 + 42393k 2 - 9161.x + S37.9 



-Poly.(l) 
-Linear [2] 



V = -1. 5Q5x + 13,35 



P = b/a 

Figure 6: Graph of SSFS, P = b/a (Note: y = k and x = Aspect Ratio, p) 

The free edge of the panel is on X - X axis. With respect to length (a), the natural frequency decreases as the 
aspect ratio increases. This model consist of one polynomial equation y = 23436x 6 - 87611x 5 + 13164x 4 - 10152x 3 + 
42393x 2 - 9161.x + 837.9 which stop at P= 1 and linear equation y = -1.505x + 13.35. 
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Figure 7: Graph of SSFS, P = a/b (Note: y = k and x = Aspect Ratio, p) 

The free edge of the panel is on Y-Y axis. With respect to width (b), the natural frequency decreases as the aspect 
ratio increases, the power equation curve is represented by y = 13.68x _1 75 
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Figure 8: Graph of SSFS, P = b/a, (Note: y = k and x = Aspect Ratio, p) 

The free edge of the panel is on Y-Y axis. With respect to width(b), the natural frequency increases as the aspect 
ratio increases, the polynomial equation curve is represented by y = 7E-05x J + 9.878x" - 0.004x + 2.469. 
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COMPARISON OF RESULTS FOR SSFS 



The table compared the results of SSFS from previous researchers with difference in percentage. In table 1, the 
present result was compared the results of SSFS from Lessia (1973) and the present result ranges from 1.44% to 12.48% 
given an average difference of 5.11% which is ok. Table 2 and 3 displayed the difference between Gorman (1982) and 
present result with aspect ratio P =a/b and P =b/a respectively. The average difference is/are 0.0167% and -7.04% 
respectively. 



Table 1: SSFS, P = b/a 



K 


D 

ph 


a 

P = b 


Present 


Lessia 
(1973) 


% 

Diff 


0.4 


10.2713 


10.1259 


1.44 


0.6 


10.7646 


10.6712 


0.88 


1.0 


12.3436 


11.6845 


5.64 


1.5 


15.4218 


13.7111 


12.48 


Table 2: SSFS, P = b/a 


K 


D 
ph 


b 

P = - 
a 


Present 


Gorman 
(1982) 


% 

Diff 


1.0 


12.3436 


12.56 


-1.75 


1.5 


10.9729 


10.90 


0.64 


2.0 


10.4933 


10.37 


1.16 


Table 3: SSFS, P = a/b 


K 


D 
ph 


a 

P = b 


Present 


Gorman 
(1982) 


% 

Diff 


1.0 


9.3259 


12.56 


-1.75 


1.5 


6.8570 


7.379 


-7.033 


2.0 


4.9369 


5.630 


-12.33 



Other researcher's worked on aspect ratio P = 1, as shown in table 4 



Table 4: SSFS, P = 1 



p = 1 


1 


2 


3 




11.6898 


11.3731 


11.684 


A. Hasan 


Baferani, A.R. Saidi, H 


Ehteshami 



2. A.R. Saidi and S.R. Atashipour 

3. S. Chakraverty 
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